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Abstract 

A destabilization in the transfer energy flux from the vacuum to radiation, for two vacuum 
decay laws relevant to the dark energy problem, is analyzed using the Landau-Lifshitz fluctua- 
tion hydrodynamic theory. Assuming thermal (or near thermal) equilibrium between the vacuum 
and radiation, at the earliest epoch of the Universe expansion, we show that the law due to 
renormalization-group running of the cosmological constant term, with parameters chosen not to 
spoil the primordial nucleosynthesis scenario, does soon drive the flux to fluctuate beyond its 
statistical average value thereby distorting the cosmic background radiation spectrum beyond ob- 
servational limits. While the law coming from the saturated holographic dark energy does not 
lead the flux to wildly fluctuate, a more realistic non-saturated form shows again such anomalous 
behavior. 
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Oftentimes, the present state of accelerated expansion of the Universe is related to some 
mysterious dark energy sector. This is linked to the longstanding cosmological-constant 
(CC) problem II], by adding to it two distinct (but connected) difficulties: (i) the "new" 
CC problem, a puzzle of why the CC is small but non-zero, and (ii) the "cosmic coinci- 
dence" /"why now?" problem js]], a puzzle concerning the current near coincidence of the 
CC energy density, /?a, with that of matter despite they scale at different rates with expan- 
sion. It seems today that a new aspect in dealing with the CC problems lies in the landscape 
of string theory [3J, though making predictions in such a theory constitutes an enormous 
challenge 4j. 

Long before the "environmental" variable CC approach of the string theory landscape 
it was noticed that, up to some extent, a traditional running of the CC can ameliorate 
the fine-tuning problems inherent to the CC by providing a viable mechanism to efficiently 
relax it from a very high value at the early Universe to its current tiny value. It was sub- 
sequently noted that some of the running CC models could successfully mimic the popular 
quintessence models as well as shed some light on the coincidence problem, thus becoming 
viable cosmological models of dark energy of the Universe. Arguably, the most appealing 
amongst them are those whose laws for the CC running can be inferred from some under- 
lying physical theory. So, some of the most attractive dark energy models involve the CC 
running laws derived from quantum theory of particle fields on the classical gravitational 
background [5], quantum gravity and gravitational holography Q. A comprehensive list 
of phenomenological laws for the CC variation under consideration well before the discovery 
of dark energy, can be found in [8J. 

Barring a time- dependent gravitational coupling or going over some scalar-tensor 



theory 



ll| , the simplest way to achieve the infrared (IR) screening of the CC is through its 



decay into matter and/or radiation. The interesting possibility in this context, put forward 
long ago, was to consider the cosmological vacuum decay into radiation as a measure of the 
temperature of the vacuum 12|. Thus, if the interactions between the two components at 
the earliest moments of the Universe expansion were fast enough to bring them in thermal 
or near thermal equilibrium, then both -the vacuum and radiation- would share a common 
temperature for, at least some time, during the expansion. Interestingly, such a scenario 
would preclude a thermal equilibrium between the vacuum and the event horizon (which 
naturally occurs in these scenarios for most of the history of the Universe), for much of 
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the time except near the Planck time. Indeed, since the temperature of the horizon as 
given by the Hubble parameter H, and the temperature of the vacuum/radiation T, scale 
differently with the expansion, they may therefore coincide (during a radiation dominated 
epoch), H ~ T, only when T ~ Mpi. 

At a macroscopic level, the transfer of energy from vacuum to radiation (and vice versa) 
is governed by the continuity equation 

Pa + Pr + 4H Pr = , (1) 



with pn being the radiation energy. Once the equilibrium (or near equilibrium) between the 
sub-systems vacuum/radiation gets established it will remain so provided the heat capacity 
of the whole system is positive-definite. Since the radiation heat capacity is necessarily 
positive, this amounts to having the heat capacity of vacuum 



C* = T[%f) (2) 



positive, where S\ represents the entropy of a variable CC. (Admittedly, there is some 
ambiguity when taking the partial derivative in Eq. ([2]) as the volume should be kept 
constant but the latter depends on T in an expanding Universe). 

Our target in the present paper is to study the fluctuations of the flux p\ entering Eq. 
(CQ) around its statistical average value, from the laws emerging from the RG-running and 
gravitational holography, taking the macroscopic criterion that C\ > as a consistency con- 
dition. To fulfill this aim we shall employ the well known Landau-Lifshitz (LL) fluctuation 



hydrodynamic theory 131 ]. which applies to equilibrium and nonequilibrium classical statis- 
I I 

tical theory [14J ]. Our particular emphasis will be on finding a scale dependence of these 
fluctuations, with a stabilization criterion that the root mean square of the fluctuations must 
never exceed the average value of the corresponding flux. 

According to Landau and Lifshitz, if the flux iji of a given thermodynamic quantity, 
which evolves in a generic dissipative process, is governed by iji = £j Yj + 5iji and the 
entropy rate obeys S = (±Yj ?/j), then the second moments of the fluctuations of the fluxes 
are given by < SijiSyj >= (T^ + Tj^SySfa — i,). The angular brackets denote statistical 
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average with respect to the reference state (i.e., Y^jT^Yj which is supposed to be steady 
or quasi-steady and constitutes the systematic part of the flux), and the fluctuations Syt 
are considered spontaneous departures from that state, thus < Syi > vanishes identically. 
The quantities and Yi stand for the phenomenological transport coefficients and the 
thermodynamic force conjugate to the flux y^, respectively. In the expression for S the 
minus sign must be taken when the product Yi iji is negative, otherwise the plus sign should 
be considered. This theory has been successfully employed to constrain models for the decay 
of the cosmological constant into radiation and/or matter [jjj] as wel 
second order nonequilibrium phase transitions in isolated black holes 

The LL theory when applied to the decay of a variable CC includes just a single flux, 
P\(t), governed by 1 



as in the analysis of 

3- 



p A = r Y + 6p A , (3) 



where the fluctuations, 5p\, coming from the decay of the vacuum are assumed to be Gaus- 
sian, random fluctuations, with uncorrelated Fourier modes due to statistical homogeneity 
and isotropy. The thermodynamic force conjugate to the flux pa follows from combining 
the entropy production rate 

Sa = YpA , (4) 

with Eq. Q. Finally, the second moment (i.e., the root mean square) of the fluctuations of 
the flux is given by 

< Sp A (ti) 8p K {tj) >= 2 T 8{U - tj) . (5) 
For the models explored below we have that 

Pa^PkH. (6) 

1 The second law of thermodynamics implies the presence of particle creation and therefore the only flux 
to be considered is p\. 
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Near t = tpi, we are in the realm of a ~ t a cosmologies (see below) thereby both terms 
on the right hand side of Eq.([6]) vary less and less faster with time. In addition, for the 
RG-running model, both p^ and H will approach a constant value at late times but well 
before the present moment (see below). Thus, the approximate steady-state regime can be 
maintained even for t — > oo, provided Ca > 0. 

Using the above sketched LL theory for the fluctuations of the fluxes, we aim to study 
their behavior for the approach to the CC problem based on the RG [5]. The RG is a 
conventional theoretical tool for investigating quantum effects and the scale dependence of 
a certain quantity. From the viewpoint of quantum theory of matter fields in curved space 
0, Q] , the renormalizability of the theory forces the vacuum action to contain the CC term 
as well as fourth derivative terms. Then, the CC term is viewed as a parameter subject to 
RG running and therefore it is expected to run with the RG scale, usually identified with 
an expansion quantity evolving smoothly enough to comply with the cosmological data. 
In such theories, therefore, even a "true" CC cannot be set to any definite fixed value 
(including zero) owing to the RG running effects. It may be surprising, however, that the 
time-dependence of the CC may be due to quantum effects from the RG, considering the 
: amiliar quadratic decoupling of heavy matter fields at low energy. 2 The reason for this result 
2o| lies in the high dimensionality (mass 4 ) of the scaling quantity p&, with the outcome that 
the more massive a field is, the more dominant the role it plays in the running -irrespective 
of the scale. Consequently, the running becomes stronger than logarithmic, thus providing a 
dynamical, efficient relaxation mechanism for the CC to go down the tiny value we observe 
today. Further, the above scenario for the CC running, with the choice for the RG scale 



p = H, may also furnish a viable cosmological model of dark energy [2ll. |22|]. The strongest 
phenomenological constraints on the RG model (within a framework where the running 
of the CC goes at the expense of the energy transfer between vacuum and dark matter) 
have been obtained recently by analyzing density perturbations for the running CC |23] and 
considering the validity of the generalized second law for the running CC scenario 



The solutions for the RG-running vacuum decay into radiation for flat space read [2jJ, [22] , 



2 Strictly speaking, the quadratic form of decoupling can be proved in a rigorous way only for higher 
derivative terms of the vacuum action, but not for the CC term itself [1q| [ . Yet, usually, the same is 
assumed for the CC term. 
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PR = Pma 4(1 v) , 



(7) 
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(9) 



where ^ 



12tt M 



t- is a dimensionless mass parameter driving the RG running. Here M 



represents an additive mass contribution of all virtual massive particles, o = ±1 depending 
on whether the highest-mass particle is a boson or a fermion, and the zero subscript denotes 
present-day value. Accordingly, \v\ ~ 10~ 2 would signal the existence of a particle with 
Planck mass (or the existence of somewhat less massive particles with large multiplicities); 
\u\ ~ 1 would indicate the existence of a particle with trans-planckian mass; \u\ ~ 10 -6 
would mean the existence of a particle with mass at the GUT-scale, whereas much smaller 
values of \u\ would imply an approximate cancellation between bosonic and fermionic degrees 
of freedom. 

For adiabatic vacuum decays one expects the Stefan law to be preserved (in the following 
we shall omit numerical factors since it is only the functional dependence that matters). 
Using Eq. (j7]) and bearing in mind that oc T 4 , we get 



-l+v 



(10) 



and since T is not expected to increase with expansion we infer that < v < 1. 

The entropy production associated to the CC decay is given by Gibb's equation with the 
chemical potential set to zero 



5a 



T 



Vpt 



(11) 



where V ~ a 3 . From Eqs. (11) and (4) we obtain 



6 



(12) 



A combination of p\ ~ a 4(1 U ^H, from Eq. ([8]), with Eq. (J3j) yields 

T~a 5u - 8 H. (13) 



We then determine the dimensionless ratio between the second moment of the fluctuations 
and the (squared) flux 

= <5pa5pa> ^ r = 

PAPA (PA) 2 



Notice that the scaling dependence of the ratio r\ is crucial. Should r\ increase with ex- 
pansion, sooner or later the fluctuations of the flux would become larger than the statistical 
average value of the flux, signaling destabilization -i.e., the flux would exhibit an erratic, 
unphysical, behavior. This would mean that there is not longer guarantee that the fluctu- 
ations of the flux preserve the equilibrium relations of the adiabatic decay of the vacuum, 
especially the Stefan law j^. Further, this would seriously upset the black-body spectrum 
of the cosmic microwave background (CMB) beyond the limits allowed by observation. The 
opposite instance [r\ decreasing with expansion), corresponds to the usual, stable condition. 

Using H ~ a -2 ^ 1 ™^ from Eq. for a <C 1, we get i] ~ a~ 5u+2 . Then, by requiring that 
rj does not increase with expansion, we obtain 

v>\. (15) 

On the other hand, for a>lwe have that H ~ constant whereby v > follows. From 
Eqs. (ED and (HD we get 

C A ~-^a 3 \ (16) 
1 — v 
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whence for the allowed range of values for v the heat capacity of the vacuum results positive- 
definite and thermal equilibrium between vacuum and radiation is to hold for ever. 

Before proceeding, it is expedient to check whether the systematic part of the flux is 
quasi-steady. It will be whenever the time scale for the vacuum to decay into radiation is 
larger than the expansion time (i.e., p\/ | p\ | > H^ 1 ). Noting that, as follows from Eq. 
([8]), the systematic part is pa = —4 p RQ v a -4 ^ 1- ^ H, it is readily seen this is guaranteed for 
v > 3/4. 

Next, we compare our LL bound, given by (1T51) . with the existing bounds on the RG 
model at any epoch. As is observationally known, at the primordial nucleosynthesis time 



the ratio Pa/pr did not exceed 0.05 [26|. Since 

Pa i 



Pr l-v 



(17) 



as follows from Eqs. (7) and (8) for a C 1, this sharply contrasts with the LL bound (15). 

The above formulae also apply in the case of a dynamical CC scenario generically dubbed 
"holographic dark energy" (HDE). Originally derived for zero-point energies clS cL bound 



on pa, the saturated form of the HDE is usually written as 27] 



Pa = ^VL" 2 , (18) 



where L denotes the size of the region (providing an IR cutoff) and c 2 is a dimensionless 
constant. This is a very important concept since for c 2 values of the order of unity, the 
HDE model also provides a very elegant solution of the "old" CC problem. Thanks to the 
relationship between the ultraviolet (UV), pa ~ A 4 , and the IR cutoff, the holographic 
information is consistently encoded in the conventional quantum field theory. The choice 



L = H 1 is clearly the most natural and simple possibility 28|, |29|, |30[ • Then, with the aid 
of Friedman's equation we can write 



c2 — • < 19 > 



where r = p R0 /p 



AO- 
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In order to connect the above formulas for the RG case to Eq. (1TB]) . we shall write p\ 
from Eq. (jSJ) in a different fashion, namely, 



p A = C + C 2 H 2 (20) 



with 



Co = PA - ^M|^ 2 , C 2 = Co + ^M 2 Pl . (21) 



The system of equations (BOD- (BID is thus equivalent to the ©-(Ej) set. The HDE law (081). 
follows from (l20l)-(l2~Tl) by setting Co = 0. Then, Eqs. ( {TBI) and ( TT9T) are readily recovered 
-modulo, the obvious identification v = c 2 . Now, as seen from Eq. (|T9|) . the bound from the 
LL theory, c 2 > |, is respected since the ratio r is tiny today, ~ 10~ 5 . Note that although 
the HDE law, Eq. (fTB|) with (fT9T) . describing the vacuum decay, does intrinsically satisfy 
the LL bound, it does disturb the big bang nucleosynthesis scenario by a wide margin. The 
main problem is that the ratio pr/ p\ stays frozen during the whole cosmic expansion so that, 
for parameters driving the accelerated expansion of the Universe at late times, a transition 
to a radiation-dominated Universe is not feasible. Thus, the model cannot be considered 
realistic. 

A more realistic class of models, which do allow transitions between the cosmological 



eras, is provided by the non-saturated HDE concept [29|, |30|, l31(- The parameterization 
is again given by Eqs. ( 1TB1) and ( TL9T) . but with c 2 promoted to a function of cosmic time 
whence the ratio between the energy densities becomes a function of time. A criterion for 
a realistic non-saturated HDE model is to saturate the holographic bound asymptotically, 
c 2 (t — > oo) = 1, while having c 2 < 1 in the radiation-dominated era. This type of 
parametrization of p\ in a non-saturation regime is particularly appealing since it reduces 
directly to (18), where again only the genuine IR cutoff shows up. It is easy to find a 
particular model that does not comply with the LL bound and the big bang nucleosynthesis 
bound simultaneously. It is possible to find such a function c 2 , which does satisfy the above 
criterion for a realistic non-saturated HDE model, so that the RG law (!20l) - (l2~T!) become 
equivalent to the law obtained from the non-saturated HDE. Indeed, the choice 
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r{t) = r ° l-ar (l-a-^)(l-a + r )^ (22) 



with e = 4a/ (1 + r ) and a = C^Hq/ p\ reproduces a non-saturated HDE model 

P, = ^c\t)L-\ (23) 



with c 2 = , equivalent to the RG-running law given by Eqs. (|20~1) - (I2T1) . As shown 
above, this model does not satisfy the LL bound without seriously affecting the big bang 
nucleosynthesis scenario. 

In summary, two of the most popular and viable dark-energy models, based on vacuum 
decaying laws, appear largely compromised by flux destabilization. Our analysis relies on 
the assumption that, at some time close to the Planck era, the vacuum was in thermal (or 
near thermal) equilibrium with radiation. We have shown that the said equilibrium would 
persist at the time of big bang nucleosynthesis, where the vacuum is restricted to a tiny 



fraction of the total energy density [26( . If, in the meantime, the fluctuations of the energy 
flux become erratic (i.e., rj > 1), then the radiation component will no longer present a 
black-body spectrum and the CMB will get seriously distorted. For a running CC scenario 
it has been shown that it is not possible to simultaneously hinder the growth of fluctuations 
(relative to the average value) and reduce the vacuum contribution at nucleosynthesis' time 
to an acceptable level. For a saturated HDE model, the flux remain under control at all 
times. However, if the nucleosynthesis bound is fulfilled, then an accelerated expansion at 
late times cannot be achieved. The non-saturated HDE model shows identical anomaly as 
the RG model. 

We may conclude by saying that either a thermal equilibrium between vacuum and radi- 
ation did never occur or the dark energy models here considered are in need of revision. 

Acknowledgments 

This work was partly supported by the Ministry of Science, Education and Sport of 
the Republic of Croatia under contract No. 098-0982887-2872, the Spanish Ministry of 



10 



Education and Science under Grant FIS 2006-12296-C02-01, and the "Direccio General de 
Recerca de Catalunya" under Grant 2005 SGR 00087. 



[1] 
[2] 

[3] 
[4] 



[5] 
[6] 
[7] 
[8] 
[9 
[10 

[11 

[12 
[13 
[14 



See, e.g., S. Nobbenhuis, |gr-qc/06090lT] and references therein. 

P. Steinhardt, in Critical Problems in Physics, edited by V.L. Fitch and D.R. Marlow (Prince- 
ton University Press, Princeton, NJ, 1997). 



See, e.g., R. Bousso, hep-th/0610211 and references therein. 



R. Bousso and I-Sheng Yang, hep-th/0703206 



I. L. Shapiro and J. Sola, Phys. Lett. B 475, 236 (2000). 
A. Bonanno and M. Reuter, Phys. Lett. B 527, 9 (2002). 

A. G. Cohen, D.B. Kaplan, and A.E. Nelson, Phys. Rev. Lett. 82, 4971 (1999). 
J. M. Overduin and F. I. Cooperstock, Phys. Rev. D 58, 043506 (1998). 
I. L. Shapiro, J. Sola and H. Stefancic, JCAP 05 (2005) 012. 

B. Guberina, R. Horvat and H. Stefancic, JCAP 05 (2005) 001. 

V. Faraoni, Cosmology in Scalar-Tensor Gravity (Kluwer, Boston, 2004); N. Banerjee and D. 
Pavon, Phys. Lett. B 647, 447 (2007). 

M. Gasperini, Phys. Lett. B 194, 347 (1987); Class. Quantum Grav. 5, 521 (1988). 
L. Landau and E.M. Lifshitz, Mecanique des Fluides (Mir, Moscou, 1971), Chap. 17. 
J.M. Rubi, in Recent Developments in Noneequilibrium Thermodynamics: Fluids and Related 
Topics, proceedings of the Meeting, Sant Feliu de Guixols, Catalonia, Spain, 1985, edited by 
J. Casas- Vazquez, D. Jou and J.M. Rubi (Lecture Notes in Physics, Vol. 253) (Springer, New 
York, 1986). 

[15] D. Pavon, Phys. Rev. D 43, 375 (1991); J. Matyjasek, Phys. Rev. D 51, 4154 (1995). 

[16] D. Pavon, and J.M. Rubi, Phys. Rev. D 37, 2052 (1988); D. Pavon, Phys. Rev. D 43, 2495 

(1991); R-K Su, R-G Cai, and P.K.N. Yu, Phys. Rev. D 50, 2932 (1994). 
[17] N. D. Birell and P. C. W. Davies, Quantum Fields in Curved Space (CUP, Cambridge, 1994). 
[18] I. L. Buchbinder, S. D. Odintsov and I. L. Shapiro, Effective Action in Quantum Gravity (IOP, 

Bristol, 1992). 

[19] E. V. Gorbar and I. L. Shapiro, JHEP 02 (2003) 021; JHEP 06 (2003) 004; JHEP 02 (2004) 
060. 



11 



[20] A. Babic, B. Guberina, R. Horvat, and H. Stefancic, Phys. Rev. D 65, 085002 (2002). 
[21] I. L. Shapiro and J. Sola, Phys. Lett. B 574, 149 (2003). 

[22] C. E. Espaha-Bonet, P. Ruiz-Lapuente, I. L. Shapiro, and J. Sola, JCAP 04 (2004) 006. 



[23] J. C. Fabris, I. L. Shapiro and J. Sola, [gr^c /0609017. 

[24] R. Horvat, Phys. Lett. B 648, 374 (2007). 

[25] J. A. S. Lima, Phys. Rev. D 54, 2571 (1996). 

[26] R. Bean, S.H. Hansen, and A. Melchiorri, Phys. Rev. D 64, 103508 (2001). 

[27] M. Li, Phys. Lett. B 603, 1 (2004). 

[28] R. Horvat, Phys. Rev. D 70 (2004) 087301. 

[29] D. Pavon and W. Zimdahl, Phys. Lett. B 628 (2005) 206. 

[30] W. Zimdahl and D. Pavon, |astro-ph/060"6555] 

[31] B. Guberina, R. Horvat and H. Nikolic, JCAP 01 (2007) 012. 



12 



